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NYO-9601 
Abstract 

This  paper  considers  the  Vlasov  equation  without  magnetic 
fields,  linearized  about  a  state  with  no  electric  fields  or 
space  and  time  variation.   Integral  representations  for  the 
Green's  functions  are  presented  for  systems  of  one,  two,  and 
three  dimensions.   For  the  one-dlmenslonal  stable  case,  two 
asymptotic  expansions  of  the  electric  field  for  large  time 
are  presented.   Similar  expansions  are  given  for  a  specific 
example  of  a  one-dlmenslonal  unstable  plasma.   In  two  and  three 
dimensions,  asymptotic  expansions  are  exhibited  for  arbitrary 
plasmas  having  cylindrical  or  spherical  symmetry. 
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On  the  Green's  Function  for  the  Linearized  Vlasov  Equation 


I .   Introduction 

One  of  the  common  models  employed  to  study  the  nature  of 
plasma  oscillations  consists  of  the  linearized  Vlasov  equation 
without  magnetic  fields  and  the  relevant  equations  from 
Maxwell's  equations.^  '      For  the  situation  in  which  the  unper- 
turbed state  about  which  linearization  is  performed  is  constant 

in  space  and  time  and  has  no  electric  f ield^  we  present  integral 

(2) 

representations  of  the  Green's  function  for  the  equations-    In 

our  treatment  the  electric  field  vector,  velocity  vector,  and 
position  vector  all  have  the  same  number  of  components  --  the 
dimension  of  the  space.   For  "stable"  plasmas  we  derive  asymp- 
totic expansions  of  the  electric  field  for  large  time.   In  the 
one-dimensional  case  we  also  obtain  asymptotic  expansions  for 
a  specific  unstable  plasma  and  comment  on  the  general  unstable 
case.   Following  the  treatment  of  Landau,  we  start  by  examining 
the  Fourier-Laplace  transforms  of  the  equations.  However,  we 
proceed  to  study  the  solutions  of  the  equations  themselves 
rather  than  their  transforms  or  the  ensuing  dispersion  relation, 

Characteristic  of  many  problems  of  mathematical  physics, 
the  form  of  the  solution  of  the  equations  depends  on  the  dimen- 
sion of  the  space.   Also,  we  find  it  necessary  to  impose 
different  conditions  and  to  use  varying  techniques  in  one,  two, 
or  three  dimensions.   Unforttmately, this  problem  presents  still 
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other  confusing  aspects.   First,  a  natural  derivation  of  the 
Integral  representation  of  the  Green's  function  does  not  seem 
feasible  except  In  special  cases.   However,  we  prove  that  the 
proposed  Green's  functions  do  satisfy  the  equations.   Secondly, 
In  certain  cases  the  solution  Is  a  distribution,  or  a  weak 
solution,  rather  than  an  ordinary  function.   Thus,  for  sim- 
plicity we  consider  the  relatively  direct  one  dimensional 
case  In  some  detail  before  proceeding  to  the  more  complicated 
two- and  three-dimensional  cases.   In  the  next  section  we  formu- 
late the  problem  generally,  and  after  a  few  comments  we  present 
some  mathematical  material  applicable  to  all  the  various  cases. 
The  third  section,  devoted  entirely  to  the  one-dlmenslonal  case, 
contains  a  derivation  of  the  solution  of  questionable  rigor 
and  a  full  proof  that  the  solution  obtained  Is  In  fact  a  solu- 
tion to  the  problem.   In  the  fourth  section  we  derive  two 
asymptotic  expansions  of  the  electric  field  for  large  time  for 
an  arbitrary,  stable,  one  dimensional  plasma.   In  the  fifth 
section  we  consider  a  specific  example  of  a  one-dlmenslonal 
plasma  that  Is  subject  to  the  two  stream  Instability  and  we 
obtain  again  asymptotic  expansions  of  the  electric  field.   We 
conclude  that  section  with  some  brief  remarks  on  unstable  one- 
dlmenslonal  plasmas.   In  the  sixth  section  we  give  the  three- 
dimensional  representations  and  then  asymptotic  expansions  for 
the  spherically  symmetric  case.   The  final  section  repeats  much 
of  the  material  of  the  sixth,  but  In  two  dimensions. 


-  5  - 


II .   Formulation  of  the  Problem 

Under  the  hypotheses  of  the  introduction  we  may  assume 
that  the  particle  distribution  function,  f(r,  v,  t),  and  the 
electric  field,   E(r,  t),   satisfy  the  equations: 


(la)  ^  (r,  V,  t)  +  V  ^+  ^  E(r,  t)-  — ^=—  =  0 


(lb)  V-E(r,  t)  =  47re  /  f(r,  v,  t)(dv) 

(Ic)  V>^E(r,  t)  =  0  . 

In  the  equations  (l),  f  (v)  is  the  unperturbed  distribution 
function  and   e   and   m  are  the  charge  and  mass  of  the  particles. 

Without  some  statement  concerning  initial  conditions  and/or 
boundary  conditions  for  the  functions   E(r,  t)  and  f(r,  v,  t), 
we  have  not  specified  a  complete  problem.   Since  we  do  not 
prove  any  existence  and  uniqueness  theorems  for  (l)  with  appro- 
priate additional  conditions,  we  can  only  propose  what  seem  to 
be  reasonable  conditions.   From  the  formulation  of  the  problem 
we  would  certainly  expect  to  give  ari  initial  condition  for 
f(r,  V,  t).   Together  (lb)  and  (ic)  imply  that  E(r,  t)  satisfies 
Poisson's  equation,  so  that  we  might  expect  to  impose  only 
regularity  conditions  at  infinity  on  E(r,  t).   Thus  we  consider 
the  Green's  function  for  the  system  (l)  to  be  functions  satis- 
fying (l)  and  the  added  conditions 


(2a)  f(r,  V,  0)  =  g(v)5(r) 


and 


(2b)  lim  E(r,  t)  -  0 

|r|  — >  00 

t  ;^  0 


Thus  the  Green's  function  represents  the  disturbance  produced 
by  a  source  of  particles  at  r  =  0   and   t  =  0  with  the  velocity 
distribution   g(v).   We  may  express  the  solution  to  the  general 
Initial  value  problem  In  terms  of  convolution  Integrals  of 
f(r,  V,  t)  or  E(r,  t)  and  the  given  Initial  data. 

While  we  have  set  up  the  equations  for  a  one  -component 
plasma,  the  system  of  equations  (l)  covers  the  many-component 
case.   For  let  us  suppose  that  there  were  many  species  of 
particles,  each  with  its  own  distribution  function,  f .  ("r,  v",  t), 
unperturbed  distribution  function   f  . (v )  ,  charge   e.,   and 
mass   m. ,   such  that  each  distribution  function  satisfied  (la). 
The  appropriate  form  for  (lb)  would  be 


(lb')  VE(?,t)  =  ^TfX  e.  /  f.  (r,  v,  t)(dv) 

(i) 


If  we  define 


f  '  (T,  V  ,  t )  =  I  ^  e^f .  (r ,  -,  t ) , 


and 
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2 
o^  ^     2  ^^  m.    oi 


then  f'(r,  v,  t),  f'(v),   and  E(r,  t)   again  satisfy  (l). 
Further,  If   f  . (v )  >  0,   then   f'(v)  >  0^  so  that  the  problem 
Is  essentially  the  same  as  with  only  one  species  of  particles. 
It  Is  well  known  and  Immediately  verified  that  the  system 
of  equations  (l)  Is  Invariant  under  the  Gallllean  transformation 

U  =  V  -  w 
r'  =  r  -  wt 
t'  =  t. 

Ultimately  we  expect  to  obtain  asymptotic  expansions  of  the 
electric  field  and  we  would  certainly  expect  them  to  vary  from 
one  coordinate  system  to  another.   For  reasons  apparent  later 
we  choose  one  particular  system  as  rest  frame,  the  unique 
system  in  which 


(5)  /  vf^(v)(dv)  =  0 


In  certain  circumstances  we  must  add  an  additional  condi- 
tion In  order  to  obtain  the  Green's  function.   The  restriction 

(^)  /  g(v)(d7)  =  0  , 

asserts  that  the  perturbation  does  not  change  the  amount  of 
charge  in  the  system..   For  example,  in  one  dimension  we  derive 


from  (lb)  that  jj—  (e(+oo,  t)-E(-oo,  t))   equals  the  charge  added 
by  the  perturbation.   Thus  if  (2b)  is  to  hold,  then  we  must 
impose  (4).   In  two  dimensions  we  do  not  require  (4)  to  hold. 
For  technical  reasons  we  start  by  imposing  (4)  in  three  dimen- 
sions, but  we  finally  drop  the  condition. 

Essentially  all  the  rest  of  this  paper  considers  each 
dimension  separately.   However,  we  present  below  certain 
suxilliary  material  which  is  independent  of  the  dimension. 
Since  we  expect  to  employ  Laplace  transform  techniques,  we  must 
study  certain  functions  analytic  in  a  half -plane.   If  we  define 
the  physical  constants 


(5a) 

and 

(5b) 


n  =  /  fQ(v)dv 


2    47rne' 

03   =  

p     m 


the  plasma  frequency. 


and  introduce  the  function   f  (v)  =  f  (v)/n  ,   then  we  define 

o  o 

the   functions 


(6) 


and 


(7) 


G(z,    k)    =    -47re 


K(v)(dv) 

/s     

k- v-z 


F(z,    k)    =  cd; 


-   (dv) 


P 


)v 


Im  z  >   0 


,        Im   z   >   0    , 


k  •  v-z 


where  k  is  a  real  unit  vector  in  the  relevant  space.   In  one 
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dimension   k   is  omitted  from  (6),  (7)  and  the  formulas  below. 
We  see  immediately  that  F(z,  k)    and  G(z,  k)  are  analytic  in  the 
domain   Im  z  >  0  if  g(v)  and  Sf  (v)/^v   are  absolutely  integrable 
Further,  if  g(v)  and  hf    (v)/Sv   have  derivatives  which  are  con- 
tinuous  and  absolutely   integrable,  then  F(z,  k)  and  G(z,  k) 
may  be  extended  continuously  up  to  the  line   Im  z  =  0.   For 
many  of  the  asymptotic  expansions  considered  below  we  must 
assume  that   f  (v)  and  g(v)  have  many  continuous  absolutely 
integrable  derivatives,  so  that  F(z,  k)  and  G(z,  k)  may  also 
possess  many  continuous  derivatives  in  the  domain  Im  z  >  0. 
However,  since  we  never  assume  analyticity  of  g(v)  or  f  (v), 
we  cannot  extend  F(z,  k)  or  G(z,  k)  analyticitally  into  the 
region  Im  z  <  0. 

We  also  need  asymptotic  expansions  of  F(z,  k)  and  G(z,  k) 
near  infinity  in  the  half-plane  Im  z  >  0.  Without  full  proof, 
from  the  identity 

r 

G(z,  k)  =  "~" 


^ire 


:(v)(d7)  1  +  X^+  -^ 


-  <^n2 

+ 


2 
z 


+ 


( V  •  k )     (v  •  k)" 


n 


z  (z-vk) 


we  see  that  if  g(v")  and  its  derivative  vanish  as  fast  as  1/ 1  v  | 

for  n   sufficiently  large,  then  we  may  derive  the  asymptotic 
expansion 

A    1    (  a, (k)       a    (k)  A  (z,  k) 

(8)          G(z,  k)=iU  +  J.l-i+...+  ^lLl4  ^11- 

^  '                    z    o     z  n        n 

I  z        z 
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valid  for  Im  z  >  0.   In  (8)^  l^n^^'  ^^  '  ^^   uniformly  bounded 


for  Im  z  >  0  and 


(9) 


a  (k)  =  ^ire 


(d7)g(v)(k.v)" 


When  g(v)  is  restricted  so  that  (4)  holdSjWe  infer  that  a^  =  0. 
In  a  similar  manner  we  obtain  the  expansion  for  F(z,  k). 


(10)      F(z,  k) 


03^       P  (k)    P^(k) 
z        z       z 


+  - 


P  (k)    B  (z,k) 
^n^  '  n^  '  ' 


n 


n 


valid  in  Im  z  >  0  and  with  |b^(z,  k) |  uniformly  bounded  there 

A 

The  coefficients  p  (k)  are  defined  by 

n 


(11)      P^(k)  =  j  {k.vf   ?^(v)(d7) 


We  have  used  (5)  to  set   P^(k)  -  0  in  (lO) 
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Ill .   The  One-Dimensional  Green's  Function 

In  this  section  we  present  a  derivation  of  the  integral 
representations  of  the  Green's  function  for  the  one  dimensional 
case.   The  derivation  is  only  appropriate  for  a  "stable" 
plasma,  but  we  show  that  the  Green's  function  obtained  is  correct 
for  any  plasma.   We  noted  above  that  we  assume  (4)  holds  so  that 
in  (8)   a   =  0.   In  this  case  we  also  ignore  equation  (ic)  which 
is  trivially  satisfied. 

If  we  perform  a  Fourier  transform  in  x  and  a  Laplace 
transform  in   t,   we  obtain  from  equations  (la),  (lb)  and  (2a) 
the  relations 

(12a)      (p  +  ikv)f(k,v,p)  +  ^  E(k,p)  ^ =  g(v),   Re  p  >  0 

and 

f 

(12b)      ik  E(k,p)  =  hire        f(k,v,p)dv  ,    Re  p  >  0. 

We  may  solve  (l2a)  for  f(k,v,p)  in  terms  of  E(k,p)  and  then 
from  (l2b)  we  get,  as  in  Landau's  paper, 

(13)  E(k,p)  =  p  G(ipA)     ^     k  >  0,   Re  p  >  0, 

k  -  F(ip/k) 

where  G(z)  and  F(z)  are  defined  by  (6)  and  (?) •   From  (13)  we 
recognize  the  well-known  dispersion  relation  to  be 

(14)  k^  -  F(ip/k)  =0,      k>0,   Rep>0. 
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According  to  the  dispersion  relation  the  plasma  Is  stable  If 
F(z)  Is  never  real  and  non-negative  for   Im  z  >  0.   Returning 
to  the  function  E(x,t),  we  derive 


(15) 


E(x,t) 


Re 


27r 


rj 


dk 


a+ioo 

-1^  pt   Ikx 
dp  e^   e 


a-ioo 


G(lp/k) 
k^-  F(lp/k) 


where  a  >  0  and  all  roots  of  (l4)  have  Re  p  <  a .   In  obtaining 
(15)  we  have  used  the  representation  for  the  delta  function 


5(x)  = 


Re 

IT 


Ikx  ,, 
e    dk 


o 


If  we  assume  the  plasma  Is  stable,  then  we  may  displace 
the  contour  of  the  p  Integral  In  (15)  to  the  left  until  a  =  0. 
We  may  now  define  the  complex  variable   z  =  Ip/k,  and  (15) 
becomes 


(16) 


E(x,t)  - 


Re 
2v-^ 


kdk 


o 


dz  ei^(^-^^)   ,Q(^) 

k  -  f(z) 


Let  us  consider 


(17) 


dz  e 


J 

—00 


Ik(x-zt)    G(z) 

k^-  P(z) 


k  <  0 


Since  the  denominator  of  the  Integrand  in  (l?)  never  vanishes 
for  Im  z  >  0,  we  may  displace  the  contour  of  Integration  in 
the  z  plane  to  any  line  Im  z  =  a'  >  0.   As  o '  —^°°,    the  Integrand 
vanishes  and  so  we  infer  1=0.   We  may  then  extend  the  k  inte- 
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^ration  in  (l6)  to  minus  Infinity  so  that 


E(x,t)  =  ^ 
27r  J 


.00  ,00 


cdk 


ds 


p 


G(^) 


ik(x-zt 


ir-  P(2) 


The   k   inteerration  Is  now  trivial  and  we  derive 


(18: 


Im 


( 


E(x,t)  =  ^  I  dz  G(z)  cos(/fTF7  (x-zt))  . 


x/t 


We  may  get  f(x,v,t)  most  easily  from  (la),  using  (2a)  and  (l8), 
and  we  find 


(19) 


f(x,v,t)  =  - 


Sf  (v) 
e   ..    o^  ^ 

—  Im  — ^^ 

mTF      dv 


dz  G(: 
x/t 


sln(/F(z)(x-zt)) 
/FTT7  (v-z) 


+  g(v)5(x-vt). 

We  obtained  (l8)  and  (l9)  only  under  the  severe  assumption 
of  stability  of  the  plasma  and  we  omitted  num.erous  troublesome 
technicalities.   Now  we  show  generally  that  (l8)  and  (19)  are 
the  Green's  function.   First  we  must  show  that  (l8)  and  (19) 
are  well-defined,  next  that  they  satisfy  the  boundary  condition 
and  initial  condition,  and  finally  that  they  obey  the  differ- 
ential equations.   We  know  F(z)     is  analytic  in  the  domain 
Im  z  >  0,  but  /f(z)    might  not  be  single  valued  in  the  domain 
or  analytic  at  all  interior  points.   Since  the  integrands  of 
(l8)  and  (19)  are  even  functions  of  )'F(z)    we  see  that  in 
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their  dependence  on  P(z)  they  are  analytic  functions  of   z   in 
Im  z  >  0.   Prom  the  asymptotic  expansion  for  F(z)  we  see  that 
for  fixed  x  and   t   the  arguments  of  the  sine  and  cosine  in 
(18)  and  (19)  are  uniformly  bounded  in  the  domain.   From  the 
asymptotic  expansion  for  G(z)^  (with  a   =  O),  we  infer  that 
the  integrals  in  (18)  and  (19)  converge  absolutely  and  uniformly 
near  infinity.   Further,  the  contours  of  integration  in  (18) 
and  (19)  may  be  taken  from  t-  to  any  point  at  infinity  in  the 
half  plane  Im  z  >  0.   We  have  just  shown  that  E(x,t)  is  well- 
defined,  but  to  show  that  f(x,v,t)  is  well-defined  we  must  show 
one  more  thing:   that  the  singularity  in  the  integrand  of  (19) 
at  z  =  V  has  no  effect.   For  y^  ^     ^^   "^ay  choose  the  path  of 
integration  so  that  it  does  not  pass  through  z  =  v  and  for 
r-  =  V  there  is  no  pole  of  the  integrand  since  the  numerator 
and  denominator  both  have  a  first  order  zero  at  z  =  v.   Thus 
f(x,v,t)  is  well-defined  and  with  little  extra  effort  we  might 
even  show  that  (f(x,v,t)  -  g(v)6(x-vt))  is  a  continuous  function 
of  its  arguments. 

We  verify  first  the  easier  of  the  initial  and  boundary 
conditions;  corresponding  to  (2a)  we  show  that 


(20)  lim   (f(x,v,t)  -  g(v)6(x-vt))  =  0 

t— ^0 


If  X  ;^  0  then  (20)  is  trivial  as  the  domain  of  integration  in 
(19)  becomes  smaller  and  smaller  as  t  -^  0;  if  x  =  0,  then  the 
integrand  tends  to  zero  uniformly  as  t  — ^  0  and  again  (20) 
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follows.   To  show  that  (2b)  holds  is  more  complicated.   We  choose 
a  path  of  Integration  outside  the  circle  |z|  =  | x/t | .   From  the 
asymptotic  expansion  for  F(:3)  in  the  form  F(z)  =  A(z)/z  ,  with 
A(z)  uniformly  bounded,  we  see  that  on  the  path  of  integration 
I /F ( z )  x|  <  ct,  for  some   c   independent  of  x  and   t.   Thus 
it  follows  that 

00 

|E(x,t) I  <  c"  j  |g(z) I |dz|  e^  ^  , 
x/t 

where   c    and   c"   are  constants  independent  of   x  and   t. 
From  the  asymptotic  expansion  for  G(z),   we  obtain  an  acceptable 
form  for  (2b) 


(21)       |E(x,t)  <  ^"^^"^ 

I  xl 


It  remains  to  show  that  the  differential  equations  are 
satisfied,  and  again  we  start  with  the  easier,  (lb).   From  (l8) 
and  (19)  we  derive  immediately  that 


(22)  Ifiii^   -   ifTre   j   f(x,v,t)dv   =    -  ^  Irn  G(|)    -   ^  g(f) 


But  from  the  integral  representation  for  G(z)  and  well-known 
formulas  we  infer  that  Im  G(x/t)  =  -^tt  e  g(x/t).   Thus  the 
right  hand  side  of  (22)  vanishes  and  (lb)  holds. 

If  we  were  to  try  to  verify  directly  the  remaining  equation. 
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(la),  we  would  find  that  f(x,v,t)  as  given  by  (19)  is  not 
differentlable  with  respect  to   x  and   t   in  the  whole  x,v,t 
space.   Since  we  have  obtained  a  Green's  function  we  might 
expect  it  to  be  a  distribution,  which  need  not  be  differentlable 
in  the  ordinary  sense.   Conveniently  we  may  show  (la)  holds 
without  any  of  the  complications  of  distribution  theory.   If 
we  define  new  variables   x'   and   t'   by 


(23a) 
(23b) 


X'  =  x-vt 

t '  =  t. 


then    (la)    becomes,    written   in   terms   of   both   the    old   and   new 
variables , 


(24) 


M^uXAl+   £  E(x,t) 
dt '  m      ^         ' 


^f    (v) 
o^    ^ 


0 


If  we  define  f'(x',v',t')  =  f(x,v,t),  then  (19)  becomes 


(25) 


f ' (x' ,v' ,t ' ) 


Sf  (v 

—  Im  3 

mw  ov 


dz 


J 


G(z)  sin(/F(z) (x'-(z-v)t ' ) ) 


v-z 


X 


7+v 


+  g(v)5(x' ) 


We  may  now  readily  verify  that  f'(x',v',t')  is  continuously 
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dlf ferentiable  with  respect  to  t'  and  that  (24)  holds.   In  a 
weakened  sense  we  have  thus  shown  that  (la)  Is  satisfied.   Hence 
(18)  and  (19)  are  the  solutions  to  the  system  of  equations  (l) 
and  {2)j   (with  (4)  assumed). 


IV.   Asymptotic  Expansions  of  the  Electric  Field  in  the  One- 
Dimensional,  Stable  Case. 

Although  we  now  have  an  integral  representation  for  the 
electric  field  in  one  dimension,  we  cannot  learn  very  much 
from  the  formula  as  it  is.   Hence,  in  this  and  the  succeeding 
section  we  obtain  asymptotic  expansions  of  the  formula  (18). 
Here  we  restrict  ourselves  to  stable  plasmas,  for  which, 
according  to  the  preceeding  discussion,  P(z)  is  never  real  and 
non-negative.   We  consider  two  separate  limits,   x   fixed  and 
t  — >  CO,   and   x/t  =  V   fixed  and  t  — >  co.   The  meaning  of  the 
first  limit  is  clear;  in  one  interpretation  of  the  second  limit 
we  move  out  on  a  possible  wave  front  traveling  with  speed   v. 
Alternately,  after  a  Gallllean  transformation  into  a  coordinate 
system  moving  with  velocity  v   with  respect  to  the  rest  frame, 
the  limit  x'  fixed,  t'  — ^  co   is  just  the  same  as  the  second  one 
above.   Because  ordinary  saddle-point  methods  for  obtaining 
asymptotic  expansions  are  not  valid  for  a  stable  plasma. 
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we  must  employ  the  Involved  procedures  given  below.   We  consider 
first  the  easier  limit,  in  which  x/t  =  v  is  held  fixed  as  t  — >  oo, 

Since  F(z)  is  never  real  and  non-negative  we  may  define 
yF{z )    as  a  single  valued  analytic  function  in  the  half-plane 
Im  z  >  0  and  we  may  also  specify  that  Im  /P ( z )  <  0  there.   If 
we  set 

(26)      j6{z)   =  j6{z,v)   =    (v-z)/F(Tr  , 

then  on  the  line  Im  z  =  0 

(27a)     Im  j6{z,v)   >  0    if  z  >  v 

and 

(27b)      Im  j6{z,v)   <  0    if   z  <  v. 


F,or  simplicity  we  temporarily  assume  that  J^       —   never 
vanishes  on  Im  z  =  0.   We  remove  this  restriction  below. 

After  these  preliminaries  we  may  now  easily  obtain  the 
desired  expansion  by  repeated  integrations  by  parts.   On 
rewriting  (l8),  we  have 

j  00 

(28)       E(x,t)  =  E(vt,t)  =  g   f  dz  G(z)e^^^(^''') 

—  00 


+  Jdz  G(z)e-"^(^'^) 


V 
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After  the  first  integration  by  parts  we  derive 


(29)    E(vt,t)  = 


Im 
27r 


Gfz)    itj25(z) 
itg'(z)  ^ 


-it^(z 
-it/)'  (z) 


+  Gl^e 


V 


V 


Im  1  J^ 
2r   1  it 


'   1  ^ 


it 


V 


We  have  obtained  (29)  under  the  tacit  assumption  that  the  two 
infinite  integrals  converge.   Prom  the  asymptotic  expansions  of 
F(z)  and  G(z)  we  see  that 


n(^)  0^2+  ttp/z  +.  .  . 

W^^   (v  +  ^2/z)[-a,^) 

Provided  v  /  0  we  obtain   (G(z)/j25  '  (z ) )  '>'   c/z    so  that  the 
integrands  in  (29)  behave  just  like  those  in  (28).   Moreover 
we  infer  from.  (27)  that  the  two  exponentials  in  (28)  or  (29) 
are  uniformly  bounded  in  absolute  value  independently  of   t. 
The  properties  just  described  hold  after  any  number  of  integra- 
tions by  parts,  so  that  after  three  more  times  we  get 


(30) 


E(vt,t) 


d_ 
dz 


;5(z,v) 


+ 


t 


z-v 


We  may  obtain  as  many  terms  in  the  expansion  as  we  please 
provided   f  (v)  and   g(v)   are  sufficiently  dif ferentiable .   We 
note  that  more  than  the  moments  of   f  (v)   and  g(v)   appear  in 
(30). 


-  19  - 


In  order  to  relax  the  restriction  that  -^   (z>v)  /=   0      on 
Im  z  =^  0,  we  remark  first  that  if  dj25/dz  vanishes  at  z  =  v 
then   F(z)  -  0,  contrary  to  the  assumption  that   F   is  never 
real  and  non-negative.   Since   d/i/dz  •^^  -  03  v/z   we  see  that 
near  infinity  d^/dz   /  0   on   Im  z  =  0.   Thus  we  may  enclose 
the  zeros  of  df)/dz      in  intervals  bounded  away  from   z  =  v 
and   z   infinite.   We  may  divide  the  integrals  in  (28)  into 
integrals  over  the  special  intervals  and  the  remainder  of  the 

line.   On  the  special  intervals  the  absolute  value  of  the 

-at 
integrand  vanishes  at  least  as  fast  as   e     as   t  — >  co  for 

some  fixed  a  >  0.   Thus  these  terms  vanish  exponentially  as 

t  — ^  00.   On  the  remainder  of  the  line^  integration  by  parts 

reproduces  the  same  series  (30). 
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The  asymptotic  expansion  for  x   fixed  and  t  — ^  oo  is  much 
more  complicated.   In  some  sense,  infinity  is  a  stationary  point 
of  the  term  in  the  exponent    z  )/'F ( z )   and  we  must  include  it 
properly.   To  study  this  lim.it  we  must  consider  the  function 
i'{z)    =   z   ^P(z)  ,  analytic  in  the  half-plane  Im  z  >  0.   From  the 
asymptotic  expansion  of  F(z),  we  obtain 


(31) 


P 


^(z)  ~  03^(1  +    ^    +    ^   + 


5 


2 


2z 


3 


Since  the  first  two  terms  of  the  asymptotic  expansion  of  ii{z) 
are  real  for  z  large  and  on  the  positive  imaginary  axis,  we 
infer  that  there  is  a  curve  on  which  i/iz)    is  real  near  the  posi- 
time  imaginary  axis  near  infinity.   This  curve  cannot  have  both 
ends  at  infinity,  for  then  the  harmonic  function  Im  ifiz)   --^ 
Im(PpCD  /2z  )   would  vanish  on  the  curve  and  at  infinity^  so 
that  Im  'i^(z)  =  0,  an  impossibility.   The  curve  can  have  its  end 
not  at  infinity  only  on  the  real  axis,  and  only  z  =  0  is 
acceptable  on  that  axis.   By  similar  arguments  there  can  be  no 
other  curves  on  which  Im  ''/'(z)  =  0  and  this  curve  cannot  cross 
itself.   We  sketch  the  situation  in  Figure  1  where  we  also 
define  two  contours,  ^-,  and  £^. 

We  start  from  essentially  the  same  formula  as  (28): 


/^^N    ^/   ,  V    Im)  I  -,  ^f    \    iKF(z)(x-zt 
(32)    E(x,t)  =  5— {    dz  G(z;e  ^     ^    ^^ 


27r 


r 


+ 


x/t 
on  i 


dz  G(z)e 


■±y¥W){x-zt 


1 


x/t 
on  £,r 
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R 


^2  \ 


Im  1/{z)    >  0 


\ 


'2 


P^-^ 


\ 


r\ 


^ 


\ 


Im  ^(z)  -  0 
/ 


/ 

/P 


7/  "\ 

rr 

Im  -y^(z)  <  0 

1 


-A 


/  ■ 


2^ 


-^    <   ' 


x/t 


Figure  1 
Curve  Im  if{z)    =   0  and  definition  of  the  contours  £,-.    and  ip 


On  Jl    ,      Im  z/ F ( z )   is  negative  so  that  the  exponential  factor 
in  the  first  integral  is  a  decreasing  function  of   t.   Similarly 
from  zero  to  infinity  on   ip   the  second  Integrand  is  a  decreas- 
ing function  of   t   while  from  zero  to  x/t  on   ^p   the  real  part 
of  the  exponent  of  the  integrand  is  negative.   Thus  any  error 
terms  in  the  asymptotic  series,  which  are  of  the  form  of  (32) 
are  all  uniformly  bounded  Independently  of   t. 

We  use  a  simple  formal  device  to  separate  the  contribu- 
tions to  (32)  from  the  neighborhood  of  infinity  and  from  the 
neighborhood  of   z  =  x/t.   We  define  an  infinitely  different labia 
function  of  one  variable  on  the  path   i-,  -  ip ,  G(z),  such  that 
G(z)  equals  the  first   n   terms  of  the  asymptotic  expansion  of 
G(z)  when   z   Is  on  £,      from  Q-.   to  infinity  or  on   i^  ^'^'^^ 


-  21  - 


Qp   to  Infinity.   Further  we  assume  that  G(s)  vanishes  Iden- 
tically from  P-.   to   x/t   to  zero  to   ?„ .   We  may  select   Q-, 
and  Qp   as  far  out  as  we  please  and   n   as  large  as  we  please, 
but  we  may  not  let   n   be  Infinite  or  Q-,   or  Qp   be  at 
Infinity.   If  we  set  t?(z)  =  G(z)  -  G(z),   then  on   i  -  I    , 
where  "3'(z)   Is  defined,  ~3-(z)  =  0(l/z'^'   )   and  also 
^(z)  =  G(z)   near   z  =  x/t. 

We  may  now  attain  the  separation  of  contributions  from 
z  =  x/t   and  infinity  by  writing 


(33a)      E(x,t)  =  E^^^(x,t)  +  E^^''(x,t) 


where 


(33b) 


E^^)(x,t 


Im 
27r 


dz  ^^^)^^y^^^^--^) 


x/t 
'on  a 


1 


+ 


dz 

x/t 
on  l^ 


^(,)e-i/PT^(>^-2t) 


and 


(33c) 


E(2)(,,t) 


g   j  dz  G(z)e^>^^^'^(^-^t) 
Vt 

-lKPlTr(x-zt)  '' 


+ 


dz  G(z)e' 


x/t 
on  I, 


n+l 


Since  ^{z)    is  0(l/z    )  we  may  Integrate  by  parts  several  times 
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In  equation  (33b)  and  we  obtain 


E(^)(x,t) 


Im 

TTt^ 


fofz)  \    1 


+  0 


(1/t^) 


z==x, 


/t 


Omitting  contributions  vanishing  exponentially  we  may 
take  the  first  (second)  Integral  In  (33c)  from  Q,   (Qp)  to 
Infinity  on  l^     (-^o)*   0'^  'the  resulting  contours  of  Integration 
In  (33c)  we  have  chosen  S'(z)  to  be  the  first   n   terms  of  the 
asymptotic  expansion  of  G(z).   Without  giving  a  full  proof  we 
wish  to  show  that  we  may  now  Insert  the  asymptotic  expansion 
for  i^(z)  In  the  exponent.   In  essence  the  procedure  Is  justified 
because  the  two  leading  terms  In  the  expansion  of  l''^(z)  produce  a 
decreasing  exponential  the  real  part  of  which  (for  Q  large  enough) 
Is   larger  In  magnitude  than  the  error  term.   On  Introducing 
the  first  few  terms  of  the  expansion  of  tA(z)  and  /F(z ) ,  we 
obtain 


(3ta) 


e(^)(x 


1    00  00 

,t)-|^|j    dz  )^(z,x,t)  +  J    dz7^z,-x,-t) 


^1 

on  i- 


on  Ir 


where 


(3^b) 


'y(z,x,t)  =  e 


■ICD  t(l  +  ^J2yi)/  a, 
D  2'      /I 


4-   2    3  , 

z     z     z 


IXO) 


2  2 

X  CD 


10) 


1  + 


o^  +  -4  (PoX-P^t)  - 


2z 


3 


3 


.  3  3 
IX  co-^ 


X 
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2 

+  -^    0,x-t(P   -  P^/4)) P  (P  X  -  P  t) 

2z    D  ^  ^  2x         ^     :? 

+  ^  [(P4  -  p'/^)-  -  (P5  -  4^  )*]  -  ^  (Pjx-t(P,-p2/4)) 

.22  2  p2^2 

IX  OJ  03   p  t 
^  Op^  -  P.t)   -    P    g     +   ... 


We  may  now  extend  the  integrals  to  zero  by  adding  parts  which 
vanish  exponentially.   The  Integrals  are  now  direct  and  we 
obtain: 


(35)      E'^^^(x,t)  =    ^    <a.cos(a3^t  +  ^ir/h) 

/^^^rp^tl  1        p 


+ 


cos(c.^t  4-  V4)  /    p^   ^2  .    ,  ^^^p"l 
o^pPgt         ^2  ^p^  -  IT  ^  ^1  +  —2 

^  \         sin(cu  t   +   7r/4)   /  3x03  6   a 


-  xcOpa2U+  0(l/t5/2)  . 


f2) 

In   E^  '^(x,t)   we  see  the  net  effect  of  Landau  damping  in 

the  reciprocal  power  decay  in  time.   The  long  range  and  collec- 
time  nature  of  plasma  oscillations  is  apparent  in  the   x 
Independence  of  the  leading  term  of  the  expansion. 
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v.   Asymptotic  Expansions  of  the  Electric  Field  for  a  One- 
Dimensional  System  Exhibiting  the  Two— Stream  Instability. 

Continuing  the  Investigations  of  the  previous  section  we 
now  examine  a  specific  example  of  an  unstable  plasma  and  obtain 
asymptotic  expansions  for  the  electric  field.   We  consider  the 
special  case  of  a  plasma  that  exhibits  the  two  stream  Instability. 
An  unperturbed  distribution  function  of  the  form 

(36)      ?^(v)  =  -I  v^e"^ 

describes  such  a  plasma.   The  largest  part  of  this  section  Is 

concerned  with  the  properties  of  the  analytic  functions   F(z), 

2  2 

z  P(z),  (v-z)  F(z),   their  derivatives  and  their  square  roots. 

After  such  a  study  the  asymptotic  expansions  are  readily  derived. 

In  examining  the  functions,  the  questions  we  seek  to  answer 

are:   (l)  Does  the  exponent  of  the  Integrand  In  the  Integral  for 

E(x,t)  have  a  saddle  point?  (2)  if  there  Is  a  saddle  point.  Is 

It  correct  to  use  ordinary  saddle  point  evaluation  of  the 

Integral?      (3)  If  there  Is  no  saddle  point  or  If  the  ordinary 

evaluation  of  saddle  point  Integrals  falls,  how  Is  the  Integral 

approximated?   We  find  that  there  is  always  a  saddle  point  but 

in  certain  cases  the  saddle  point  contribution  is  negligible 

compared  with  the  end  point  terms. 

From  the  definition  of  F(z)  and  (36)  we  obtain  easily  that 


(37)      F(z)  = 
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If  we  set   z  ==  X  +  iy,  then 


(38) 


-u 


u-z 


du 


^^  -  X  +  ly -u 

(u-x)  +  y 


From  (37)  we  see  that  on  the  real  axis  Im  F(z)  =  0  only  at 

z  =  t  1  and  at  zero,  and  at  these  points  Re  F(z)  is  negative 

and  positive  respectively.   Equations  (38)  and  (57)  imply  that 

F(z)  Is  real  on  the  Imaginary  axis.   Since  F(z)  is  real  and 

positive  at  z  =  0  and  real  and  negative  on  the  imaginary  axis 

near  infinity,  it  follows  that  at  some  point   z    on  the 

o 

Imaginary  axis  F(z  )  =  0.   Further,  since  F(z  )  =  0  and  at 
infinity  F(z)  =  0,  Rolle's  theorem  implies  that  F'(z,)  =  0  for 
some  Z-,  on  the  imaginary  axis  between  z   and  infinity.   We  use 
the  principle  of  the  argument  to  count  the  zeros  of  F(z)  In  the 
upper  half  plane. 

As  the  variable  z  follows  the  contour  of  Figure  2a,  the 
argument  of  F(z),  shown  in  Figure  2b,  changes  by  27r.   Hence, 


Figure  2a 

Definition  of  the  contour  used 
to  count  the  zeros  of  F(z) 


Figure  2b 

Path  of  the  argument  of  F(z)  as 
z  follows  the  contour  of  Fig. 2a 
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we  infer  that  F(z)  has  one  zero  In  the  upper  half  plane  and  it 
must  be  at  z  . 

Prom  the  equation  for  F(z)  we  derive  immediately  that 


09) 


F'  (z) 


fco 


CO  2 


p 


2    ^    4x         e' 


j/? 


2yT  (z-'-z)  +  (-1+  5z    -2z^) 


u-z 


du 


On  the  real  axis  Im  F'(z)  vanishes  at  the  four  real  roots  of 

4    2 
2z  -  5z  +1=0.   For  z  real  and  positive  Re  F'(z)  is  positive 

at  the  larger  root  of  the  polynomial  and  negative  at  the  smaller 

root.   For  z  real  and  negative  the  same  relations  hold.   To 

count  the  zeros  of  F'(z)  we  use  Figures  3a  and  3b. 


Figure  3a. 

Definition  of  the  contours  used 
to  count  the  zeros  of  F'(z) 


Figure  3b 

Path  of  the  argument  of  F'(z) 
as  z  follows  the  contour  of  Fig. 3a 


Thus  F'(z)  has  one  zero  in  the  upper  half  plane  and  z-,  is  that 

F(z) 


zero . 


It  is  useful  to  examine  whether  Im  ^ )  7 /.  \  ever  vanishes  on 
the  real  axis.   If   Im  Z!)y  v  =  0  then  Re  F(z)  Im  F'(z)  - 
Re  F'(z)  Im  F(z)  =  0.   From  equations  (37)  and  (39)  we  infer  that 
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Im  -p)(\    =  0  on  the  real  axis  only  if   z  +  1  =  0.   Hence, 

Im  „  ;7\   is  of  one  sign  there. 

We  are  now  prepared  to  draw  several  contour  plots  of  F(z) 

and  related  functions.   In  Figure  4  we  plot  the  curves  on 

which  Im  F(z)  =  0. 

I 
I 
I 

F(z)  real  and 

positive 

F(z)  real  and 
negative 


/ 

/ 
/ 

1 

1 

1^ 

1 

Ir 

-1 

0                 +1 

Figure   4 

Curves   Im  F(z)    =   0 

We  know  that  F(z)  is  real  on  the  imaginary  axis  and  it  is 
positive  or  negative  there  as  shown  in  Figure  4.   Since 
F'(z  )  =  0,  four  lines  leave  z   on  which  F(z)  is  real.   By 
symmetry,  the  two  lines  not  yet  accounted  for  either  both  go 
to  the  real  axis  or  to  infinity.   Since  Im  F(z)  is  harmonic 
and  vanishes  at  infinity,  the  two  lines  cannot  go  to  infinity 
of  Im  F(z)  =  0.   Thus  the  lines  leaving  z-,  connect  as  shown  in 
Figure  4.   By  similar  arguments  there  are  no  other  curves  on 

which  Im  F(z)  =  0. 

2       2 
Next  we  consider  the  function  if    (z)  =  z  F(z).   On  the 

2  2        2 

imaginary  axis  if    (z)  is  real  and  ii    (o)  =   if    (z  )  =  0.   For  some 

z„  on  the  imaginary  axis  between  z   and  zero  -r—  if    (z)  =  0. 
2  D    J  Q  dz    ^  ' 
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Just  as  for  F(z)  we  obtain  the  contour  plot  of  Figure  5' 


V'  (z)  real  and 
positive 

if    (z)  real  and 
negative 


i^l 


*  o 


,--"1*2. 


-1 


+  1 


Figure  5 
2 
Curves  Im  V'  (z)  =  0 


In  treating  the  stable  plasma  we  introduced  the  function 
)6{z,y)    by  a  relation  equivalent  to 


(40) 


J252(z,v)  ^    (v-zfFiz) 


We  are  again  interested  in  this  function,  and  in  order  to  use 

saddle  point  methods  we  must  find  values  of   z   for  which 

d    2 

-^  j^    {z,v)    =   0.   Such  values  satisfy  the  equation 


2F(z)  , 


We  are  naturally  lead  to  introduce  the  function 


(41) 


X(z)  =z.|^  =  ^^f^^2(^) 
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such  that  the  real  values  of  '^(z)  correspond  to  points  z'  and 
values  of  V  =  ^(z')  such  that  ^  j^    (^Sv)  =  0.   On  the 
imaginary  axis,  where  A.i^)    is  defined  It  Is  Imaginary.   Thus 
curves  on  which  Im  y^{z)    =   0  can  cross  the  Imaginary  axis  only 
at  Z-,  (where  'X(-)  has  a  simple  pole)  or  Zp  (where  ^(z)  has  a 
simple  zero).   As  'V(z)  Is  not  real  on  the  real  axis,  we  may 
draw  the  curves  on  which  ImX(z)  =  0  as  In  Figure  6.   On  the 


Im  y(: 


=  0 


Figure  6 
Curves  Im')^(z)  -   0 

contour  on  which  Im'^(z)  =  0,  '\/(z)  assumes  all  real  values  at 
least  once.   We  have  not  really  shown  that  the  curves  leaving 
Z-,  and  Zp  do  not  meet  at  infinity.   Using  the  asymptotic  expan- 
sion for  ~)C{z)    we  might  actually  show  this  is  not  the  case,  but 
as  we  do  not  need  this  precise  statement  we  pass  over  the 
question . 

The  function  we  wish  to  consider  finally  is  the  one  occurring 
in  the  exponent  of  the  integral  representation  for  E(x,t),  and 
this  function  is  yj6    (z,v)  =  (v-z)  )/'f  (  z  ) .   HoweverjWe  cannot 
easily  consider  (v-z)  )/F(z )  as  it  has  branch  points  in  the  half 
plane  Im  z  >  0.   Thus  we  examine  ^    (->v)  itself  and  then  we  shall 
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be  able  to  make  any  necessary  analytic  continuations  of  ff)    (z,v). 
The  Information  we  wish  to  obtain  is:   where  Is  f)    (z,v)  real,  and 
is  it  positive  or  negative  there?   In  the  following  discussion 
we  assume  v  /  0,  +  1.   By  previous  arguments  at  least  one  end 
of  any  contour  on  which  Im  j2^  (z^v)  =  0  must  lie  on  the  real  axis. 
On  the  real  axis  ^    (z,v)  is  real  at  t  1,  where  it  is  negative; 
at  z  =  0,  where  it  is  positive;  and  at  z  =  v  where  ^    (z,v)  has 
a  double  zero.   Prom  the  points  z  =  0  and  z  ==  +  1  there  should 

be  contours  going  into  the  upper  half  plane  on  which  f)    (z,v)  is 

2 
real.   Leaving  z  =  v  should  be  two  contours  on  which  f)    (z,v) 

2  2 

is  real  and  on  one  j?5  (z,v)  is  positive,  on  the  other  ^    (z,v)  is 

negative.   Of  these  two  contours  at  most  one  can  go  to  infinity. 

The  remaining  information  that  we  have  is  that  on  the  imaginary 

p  2 

axis  f)    (z,v)  is  not  real  except  at  z  =  z   where  P    (z,v)  has  a 

2 
simple  zero.   Thus  only  one  curve  on  which  Im  j?5  (z,v)  -   0  can 

cross  the  imaginary  axis  and  the  crossing  is  at  z  =  z  .   We 

must  now  determine  how  to  connect  all  the  various  curves.   We 

make  the  connections  for  the  cases  of  v  small  and  v  large  and 

then  conjecture  a  continuous  transition  between  the  limits.   We 

use  the  properties  that  for  v  large  ^    (z,v)  is  approximately 

p  p 

V  P(z)  and  for  v  small  f)    (z,v)  is  approximately  -^{z). 

2 
In  Figure  7a  we  draw  the  curves  1v\  f)    (z,v)  =  0  for  v  very 

large;  in  Figure  Th  we  draw  the  same  figure  showing  how  the 

2 
curves  approximate  those  of  Im  v  F(z)  =  0.   To  justify  Figure 

2 
7  we  note  first  that  from  (37)  and  (39)  we  may  expand  f)    (z,v) 

near  z  =  0,  and  we  obtain 
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f{z,v)    =   v^|F(0)|  +  v^z(l|F'(0)|  -  ^^Mi)  +  0(z2 


0    {z,y)    real  and 
positive 

j6    (z,v)  real  and 
negative 
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Figure  7a 
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Figure  7b 
2 
Curves  Im  f)    (z,v)  =  0  for  v  large 


Thus,  for  V  >  0  the  curve  on  which  'jb    (z,v)  is  real  and  positive 
leaves  the  origin  and  goes  into  the  first  quadrant,  as  shown. 
Next,  let  us  consider  the  values  of  ^    (z,v)  on  a  semicircle  of 
radius  larger  than  1  and  | z, | ,  but  much  less  than  v.   If,  for 
V  large,  p    (z>v)  is  never  real  and  positive  on  that  semicircle, 
then  the  contour  line  leaving  the  origin  must  go  to  z   and  then 
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connect  to  -1.   Since  both  contour  lines  from   v   cannot 
go  to  Infinity,  they  must  connect  as  shown.    It  remains     :: 
to  show  that  on  our  given  semicircle  f)    (z,v)  is  not  real  and 
positive.   Let  us  suppose  the  contrary,  that  there  be  infinitely 


many  pairs  v.,  z.,    with  v. 


>,  such  that  f)    (z.,v.)  is  real  and 


positive.   Since  the  semicircle  is  compact,  we  may  pick  a  subset 

of  z ..  '  s  which  converge  to  some  point  ^.     As   (1  -  (rp-.'  )f(z.) 
J  ■  j'      J 

is  real  and  positive  for  all  j,  when  j  —^  0°,  z  — •  C,   and  v  .  — 7  <», 

J  J 

so  that  F(0  is  real  and  positive,  which  is  false.   Again  we 

have  not  shown  that  the  contours  from  -1  and  z   do  not  meet  at 

o 

infinity.   While  we  might  show  this  with  a  little  effort,  since 
the  result  is  not  necessary  below  we  ignore  it. 

In  Figures  8a  and  8b  we  give  similar  diagrams  for  0  <  v  «  1 


^    (z,v)  real  and 

positive 
0 
^    (z,v)  real  and 

negative 


■1         V     +1 
Figure  8a 
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Figure  8b 
2 
Curves  Im  /!>  (z,v)  =  0  for  v  small 
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To  justify  Figures  8  let  us  consider  first  a  semicircle  with 
center  at  the  origin  and  passing  thro^igh  Zp .   By  arguments  slml- 
lar  to  the  one  above.  If  f)    (z,v)  Is  real  and  positive  on  the  semi- 
circle as  V  tends  to  zero,  then  ii{'z)    Is  real  and  positive  there  -■ 
a  contradiction.   Thus  the  curves  from  the  origin  and  z  =  v 

connect.   By  taking  a  similar  semicircle  with  radius  larger 

2 
than  1  and  | z. | ,  we  may  show  that  for  v  small  P    (s^v)  Is  never 

real  and  negative  there.   Thus  the  remaining  contours  are 

correctly  connected. 

In  order  to  make  the  transition  between  v  small  and  v 

large  we  use  the  properties  of  'X(z)  given  In  Figure  6.   We 

recall  that  values  of  z  on  Im  "XC^)  =  0  correspond  to  values  of 

V  =  X(-)  such  that  -T:r  fi    (-^v)  =  0.   In  other  words,  each  point 

on  the  curve   Im  "X(z)  =  0  Is  a  double  point  of  the  contour 

2 
lines  of  P    (zjv)  for  the  appropriate  value  of  v.   From  the 

general  configuration  of  Figures  7  and  8  we  see  that  only  one 

2 
curve  on  which  f)    (z,v)  Is  real  and  negative  appears  In  either 

2 
quadrant;  thus  If  any  curves  1v\  f)    (z,v)  =  0  have  a  double  point, 

2 
then  we  expect  that  f)    (z,v)  Is  positive  there.   From  figures  6, 

7,  and  8  we  see  that  for  v  large  the  value  of  z  for  which 

^i —  ti    (z,v)  =  0  Is  near  z,  while  for  v  small  It  Is  near  z„ .   For 
dz  '^  ^    ^  1  2 

no  value  of  v  should  the  double  point  lie  Inside  the  curve  from 

2 
z  =  +1  to  z  =  V.   As  this  curve  corresponds  to  f)    (z,v)  real  and 

negative  we  expect  no  double  points  on  it,  so  that  as  v  — >  1  we 
find  that  z  =  v  =  +1  is  a  double  point  of  (v-z  )j/F(z)  --  a  con- 
tradiction.  Thus  for  V  large  the  double  point  may  be  connected 
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to  the  negative  real  axis  by  a  curve  not  passing  through  any 
curves  on  which  Im  /)  (z,v)  =  0.   For  small  v  no  such  connection 
Is  possible.   Since  all  curves  move  continuously  with  respect 
to  V  we  expect  to  find  some  value  of  v   for  which  the  double 
point  crosses  a  curve  Im  /)  (z,v  )  =  0.   Figure  9  Illustrates 

this  critical  situation.   In  Figure  9  we  have  omitted  the  curve 

o 

from  z  =  +1  to  z  =  V  on  which  fi    (z,v)  Is  real  and  negative. 


^    (z,v)  real  and 

positive 

^    (z,v)  real  and 
negative 


Figure  9 
2 
Curves  Im  ^    (z>v)  =  0  for  v  = 


V 


For  V  <  V   we  expect  the  stationary  point  of  ^    (z,v)  to  lie 
within  the  curve  from  -1  to  z   to  Infinity.   For  v  >  v  we  expect 
the  double  point  to  be  outside  the  curve  from  +1  to  v  to  Infinity. 
The  value  of  v   is  given  by  the  simultaneous  solution  of  the 
relations 


(42a) 
(42b) 
(42c) 


Im  X{^^)    =   0 


V   =  ^(z,) 
c       5 


Im  /)    (z^,v^)  =  0 
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We  show  that  the  asymptotic  expansions  change  character  as  v 

becomes  larger  than  v  .   It  is  possible  that  several  solutions 
^         c 

of  the  equations  (^2)  exist,  but  we  assume  this  does  not  happen. 
The  methods  below  would  apply  in  any  case. 

We  can  now  proceed  to  the  principal  task  of  obtaining  asymp- 
totic expansions  of  the  electric  field.   Again  we  consider  first 


the  limit  x/t  =  v  fixed  and  t 


We  define  a  special  branch 


■J 

of  the  square  root  such  that  Im  jZ!i  (z  ,v)  >  0.   We  call  f){z,\i)    the 
function  ^    (zjv)  so  defined.   If  we  consider  the  analytic  func- 
tion ^(z,v)  and  perform  an  analytic  continuation,  then  as  we 
cross  a  curve  on  which  f)    (z,v)  is  real  and  positive  j^(z,v)  moves 
onto  the  other  sheet  of  the  Riemann  surface  on  which  Im  /3(z,v)  <0, 
For  the  case  v  >  v  we  define  contours  of  integration  in 

Figure  10.   Again  we  have  suppressed  the  curve  from  z  =  v  to 

2 
z  =  +1  on  which  p    (-jv)  is  real  and  negative. 


p    (-jv)  real  and 

positive 
2 
P    (z,v)  real  and 

negative 


-1 


v 


Figure  10 


Definition  of  contours  i-,  and  ^p  when  v  >  v 
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starting  from  (l8)  for  E(vt,t)  with  i,  as  the  contour  of  Inte- 
gration, we  may  make  an  analytic  continuation  of  /){z,\r)    across 
the  curve  on  which  j6    (z,v)  Is  real  and  positive,  and  then  we 
distort  the  contour  of  Integration  to  l^   for  one  of  the  terms 
and  we  find 


(^3) 


E(vt,t)  = 


Im 
27r 


lt/!i(z,v)^/  N-i 
?  ^ "-    '     ^G(z)dz  + 


e-^^^("'^^G(z)dz 


v 
on  £- 


V 

on  £, 


In  equation  (43)  both  terms  are  decreasing  exponentials  and  we 
may  derive  equation  (30) . 

Next  we  consider  0  <  v  =  x/t  <  v  ,   v  fixed  and  t  — >  oo. 
In  Figure  11  we  define  the  relevant  contours,  again  omitting  the 
curve  from  v  to  +1. 


j6    (z,v)  real  and 
positive 

j6~{z,v)    real  and 
negative 


Figure  11 

Definition  of  contours  £^  ,    Ir.,    and  1-,   when  v  <  v 

1    d       j5  c 


We  start  from  the  same  equation  for  E(vt,t)  with  the  contour 
of  Integration  i-,  .   On  making  the  analytic  continuation  across 
the  contour  from  z   to  Infinity  we  may  write 
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(44) 


E(vt,t)  =  g|  I  G(z)e^^^("'^)d 


z  + 


e   ^^         Mz  G(: 


V 

on  a. 


o 
on  ^ 


3 


+ 


Im 
27r 


e-"^("'^)G(z)d. 


J 

V 

on  i 


2 


In  equation  (44)  the  first  two  integrals  are  decreasing  functions 
of  time  and  are  bounded  by  some  constant  independent  of  x  and  t. 
However,  the  third  integral  is  an  increasing  exponential.   Thus 
we  assert 


z 


(^5) 


E(v,t) 


Im 
27r 


e-"^^"'^)G(z)dz  +  0(1). 


J 

V 

on  i. 


From  the  earlier  remarks  we  expect  to  find  a  value  Zu 

within  the  contour  from  -1  to  z   to  infinity  at  which 

o  ^ 


d^ 


d 


-  (z,v) 


Z  =  Zi 


0.   Since  j25(z,v)  =  0  at  z  =  z   and  z  =  v,  and 


since  Im  j6{z,y)    >  0  within  the  chosen  contour  we  should  be  able 
to  deform  the  contour  of  integration  £^    so  that  it  passes  through 
z^i  and  so  that  for  large  t  the  absolute  value  of  the  integrand 
is  largest  on  ip  at  Zu.      We  may  then  perform  a  saddle  point 
asymptotic  expansion  and  we  get 

„„  (       -ityj^(z„v)  , ::, 

(46) 


Re 


E(t,vt)-^  2^<JG(z^)e  /t^"(z^,v) 


f 


IT 
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The  square  root  in  the  exponent  of  (46)  has  negative  imaginary 
part,  the  other  square  root  is  chosen  with  Re  /jO"  (Z||  ,v)  >  0. 
If  jt)^^  {zu)    is  real  and  negative  then  (46)  is  a   iguous  up  to  a 
sign  which  depends  sensitively  on  the  contour  H^.      As  v  — ^^^ 
from  below,  z^,  approaches  z-,   and  the  growth  exponent  goes  to 
zero. 

Equation  (46)  includes  the  case  x  fixed  and  t  — >  °o,  with 
the  replacements  Z|^  — >  z^,    v  — >  :^^  and  ^"{zi^,v)   — >  ii"  {z    ) . 

Qualitatively  we  see  that  the  electric  field  is 
oscillating  and  exponentially  growing  behind  a  wavefront  moving 
with  speed  v  =  v  .   Ahead  of  the  front  the  disturbance  disappears 
as  t  — >  CO,   At  a  fixed  point  behind  the  front,  as  t  — >  <»  the 
electric  field  only  grows  and  does  not  oscillate.   The  maximum 
growth  rate  exponent  behind  the  front  is  given  by  -iZo  /^(^p^  ' 

We  can  say  very  little  concerning  a  general  unstable 
plasma  as  it  is  easy  to  produce  various  examples  of  plasmas  for 
which  the  asymptotic  expansion  of  the  electric  field  has  very 
Irregular  behavior.   Perhaps  the  strongest  statement  one  might 
make  is  that  a  plasma  unstable  in  the  sense  that  the  dispersion 
relation  has  roots  for  Re  p  >  0,  will  almost  assuredly  be 
unstable  in  the  sense  that  in  some  coordinate  frame  the  electric 
field  will  grow  exponentially  in  time.   However,  the  coordinate 
frame  need  not  be  the  rest  frame  (relative  to  the  definition  (4)) 
In  fact,  the  frame  may  be  moving  with  arbitrarily  large  speed 
relative  to  the  rest  frame.   Also,  the  exponent  of  growth  will 
usually  be  less  than  that  predicted  by  the  dispersion  relation. 
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In  the  example  considered  above  the  region  of  instability  was 
limited  by  the  condition  |x|/t  <  v  .   Again  systems  exist  for 
which  there  is  no  such  v  .   With  little  effort  the  reader  may 
verify  that  for  appropriate  values  of  the  parameter  5  the 
unperturbed  distribution  functions  below  have  the  properties 
just  described. 

r^  n  (         ^  ('f;  \  \2  -(u-b(5)  ) 


where 


(47b) 


5  . 


.^.i 


(47c)     a  =  6  +  b 
and 

(48)      'f^(u)  =  4  ^-^  e-"  cos2(5x) 

^  (1  +  e-^  ) 
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VI .   Three-Dlmenslonal  Green's  Function 

Having  concluded  our  treatment  of  the  one-dimensional  case 
we  now  consider  the  three-dimensional  case.   We  obtain  the 
Integral  representations  of  the  Green's  function  by  an  Indirect 
procedure.   Instead  of  determining  functions  satisfying  the 
Initial  condition  (2a),  we  exhibit  a  Green's  function  that 
satisfies  the  Initial  condition 


(2a')      f(r,v,0)  =  ^  g(v)  , 


where   A   Is  some  constant.   In  some  sense  we  should  be  able 

to  obtain  the  desired  Green's  function  by  applying  the  operator 

1    2 
-  T, — r  to  that  solution.   In  order  to  avoid  the  addition  of  an 

47rA 

Infinite  amount  of  charge  to  the  system,  when  we  use  (2a')  we 
must  Impose  (4).   In  the  final  Green's  function  we  do  not 
require  (4 ) . 

It  has  been  conjectured  that  If  the  equation  (la)  Included 
magnetic  effects  and  If  all  of  Maxwell's  equations  were  added, 
then  the  only  stable  unperturbed  distribution  functions  would 
be  spherically  symmetric  In  velocity  space.   Since  we  have  not 
Included  magnetic  forces  we  might  expect  our  result  to  be 
relevant  to  real  physical  problems  only  when  f  (v)  Is  a  function 
of  |v| .   Thus,  after  obtaining  the  electric  field,  we  consider 
the  case  In  which  f  (v)  and  g(v)  depend  on  |v|  only,  and  we 
obtain  asymptotic  expansions  as  t  — >  oo. 
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We  assert  that  the  solution  of  the  system  of  equations  (l), 
(2a'),  (2b),  and  (^)  Is 


(49)     E(r,t)  =  Im 

J 


and 


dfi  k   dz  G(z,k)cos(yF(^,k) (k-r-st) ) 

.A.    

k  •  r 


(50)     f(r,v,t)  =  -  I  Im 


m 


av 


r  dz  G(z,k)sln(/F(z,k)'(k-r-zt)  )  [   27r^g( 
(k'v-z)yF(z,k) 


r-vt  I 


In  (49)  and  (50)  k  Is  a  three  dimensional  unit  vector  and  dfi  Is 
the  surface  element  on  t-he  k  unit  sphere.  Formulas  (49)  and  (50) 
are  natural  extensions  of  the  one  dimensional  result  and  they 
may  be  derived  from  the  transforms  as  In  Section  III.   If  we 
consider  the  integrand  of  the  k  Integration  In  (4-9)  or  (50)  we 
see  that  the  expression  Is  the  same  formula  as  (18)  or  (19)- 
Thus  the  arguments  of  the  one  dimensional  case  show  that  (49) 
and  (50)  are  well  defined.   Here  we  must  use  the  condition  (4) 
to  guarantee  that  the  Integrals  converge  at  Infinity 
has  no  other  function. 

¥e  next  seek  to  show  that  f(r,v",t)  satisfies  the  Initial 
condition  (2a')  and  E(r,t)  the  boundary  condition  (2b).   If  we 
define 
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(51) 


L(k,r,v,t)  =   dz  G(z,k) 

x^  _ 

k-  r 


((kT-zt)/F(z,k)') 


(k-7  -  z)yF(7:iT 


then 


f(F,7,t)  4-  ^fs(^) 

I r-vt I 


Sf  (v) 

—  Im  — 

m 


r 


hv 


an   k  L(k,r,v,t) 


For  r  =  0  then   11m  L(k,r,vt)  =  0,  so  that  the  Initial  condition 

t  ^>  0_ 
Is  satisfied.   For  |r|  ^   0,    we  divide  the  unit  sphere  Into  two 

regions,  such  that  In  the  first  | (k-r)|  y   t  ^     ,    and  In  the 

I  ,/v  _  ,  I     1/2 
second  | (k-r)|  <  t  ^  .   In  the  first  region  we  obtain  easily 

I  ,^  —  —      ,  I      1/2  1  ,^ \  I 

|L(k,r,v,tj|  <  ct  '  ,  and  In  the  second  region  |L(k,r,v,t)|  <  c'. 

Since  the  total  solid  angle  contained  in  the  second  region 

vanishes  as  t  — ^  0,    we  derive 

lim  (f(F,7,t)  +  2z!k(vI)  =  0. 
t->  0  I  r-vt  I 

In  order  to  derive  the  boundary  conditions  for  E(r,t)  we 

must  modify  the  asymptotic  expansions  for  G(z,k)  and  F(z,k). 

Since 

,   ^ ,      ,2 
G(z  ,k)  -v^  a-j/z  , 

we  can  find  positive  constants  a'  and  c  such  that 


(52) 


G(z,k)|  < 


a' 


I  2^ 
z   +  c 
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In  the  whole  half -plane  Im  z  >  0.   Similarly  we  have 


(53)      |F(c,k) 


P 


z  +  c 


Thus  if  Ik-r/t|  <  |z|  then  | /F (z ,k) (k • r-zt ) |  <  c"t  for  some 
constant  c".   Hence 


E(r,t)|  < 


an 


dce(+'="t). 


a 


J 

k«  r 


,2 
z   +  c 


After  some  elementary  integrals  we  derive 


(5^) 


E(r,t)|  <  c'"  e^"^  ^  |log(t/r; 


so  that 


(55; 


lim   |E(r,t) |  =  0 
I  r|  — >  00 
t>0, fixed 


The  proof  that  the  first  differential  equation  holds  is  the 
same  as  before.   Again  (la)  does  not  hold  as  it  stands,  but 
after  the  change  of  variables 


r'  =  r-vt 
t  '  =  t  , 


(la)  becomes 


(la" 


(r.v,t)  +  -  E(r,t) 


^ 


=  0, 


a^ 
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Equation  (la")  Is  essentially  trivial  to  verify.   Equation  (ic) 
Is  also  trivially  checked.   We  now  consider  (lb). 
We  derive  readily 


V»E(r,t)  -  47re    f(r,v,t)dv  =  -  ^  Im   dQ   G(z,k) 


z  =■ 


k-  r 


+  Qw^e    I  dv 


1 


:(v) 


r-vt 


so  that  we  must  show 


(56)  Im        dfi  G(z,k) 


J 


z=k>  r/t 


=   tiiPe 


d7    _gl^) 
|v-r/t) 


We  can  rewrite  the  defining  equation  for  G(z,kj  as 


G(z,k)  =  hire        dv(k-  ^  g(v  )  )log(k- v"-z) 


Hence  we   derive 


J 


-Im        dfi  G(z,k) 


z=k  •  r/t 


=    -2-/e 


dn(dv)(k.    —  g(v))sgn(k- (v-r/t)) 

av 


The  k  Integration  Is  elementary  and  (56)  follows  after  an  Inte- 
gration by  parts. 

While  we  have  obtained  the  Green's  function  for  the  equations 
with  the  Initial  condition  (2a'),  we  have  not  yet  solved  the 
original  problem  with  Initial  condition  (2a) .   We  assert  that 
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(57a)     E(r,t)  =  ^  V^^E(r,t 


and 


(57b)     f(r,v,t)  =  ^  7^"  f(r,v,t 


2_    r72  ■;i, 


^  _  _ 

are  the  desired  solutions.   Since  the  function  f(r,v,t)  is 

—    2  '^ 

not  twice  dif ferentiable  with  respect  to  r,  V  f  must  be  a 

distribution.^  ''   To  see  the  , meaning  of  a  distribution  as  used 
here  we  consider  functions  t{v)    which  are  two  times  continuously 
differentiable  and  which  vanish  at  infinity  at  least  as  fast  as 
l/r  ,  and  whose  first  two  derivatives  also  vanish  that  fast.   To 
define  the  distribution  f(r,v,t)  we  must  define  all  integrals 
/  j?5(r)f  (r,v  ,t  )dr,   but  we  have  naturally  that  /  (dr  )/i(r  )f  (r  ,v,t ) 
/V  j?5(r)f  (r,v  ,t  )dr ,   Thus  for  an  arbitrary  /i(r)  of  the 


\^ 


class  considered,  the  functions 

(£(F,t)  =  ^  /  //'(r';  E(F-F',t)(dF') 

3(r,v,t)  =  ^/V/'(r')  f(r-rSv,t)(dr' ) 


577-^ 


satisfy  the  differential  equations  and  c(r,t)  vanishes,  |r|  ^  0°, 

and   ■/(r,v,0)  =  -  ip-  /V'^.J^fr')   _  ^  _     (dr  '  )  =  /  /5(r  '  )  6  (r-r  '  )  (dr  '  )  . 

I r ' -r i  _  _ 

Hence,  in  the  sense  of  distributions,  or  weak  solutions,  E(r,t) 

and  f(r,v,t)  are  the  desired  solutions. 

The  electric  field  in  the  Green's  function  is  actually  a 

proper  function,  for  we  derive  easily 
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E(r,t) 


Im 


(2Tr) 


5 


dn  k 


dz  G(z,k)F(z,k)cos(/F(z,k) (k-r-zt)) 


■J 

k-  r 


Im  a 


r 


Q-w\  ^      u 


r  <^  ^  /k  •  r 


dfi  kkG(\^  ,  k)  , 


and  after  some  elementary  Integrals  we  obtain 


(58) 


E(r,t)  =  -^ 


(dv)g(7)  ^ 


v-r/t 


Im 


(27r)- 


dfi  k 


dz  G ( z, k )P ( z, k ) cos (yF(z,k)' (k-r-zt)) 


k'  r 


We  might  obtain  an  almost  comparable  formula  for  the  distribu- 
tion function,  but  it  is  very  cumbersome  and  of  little  interest. 
However  we  note  that  the  contour  integral  in  (58)  now  converges 
at  infinity  even  though  we  no  longer  require  that  (^)  hold.   We 
might  easily  show  that  the  same  is  true  of  the  distribution 
function,  for  if  in  (50)  we  take  the  path  of  integration  from 
k-F/t  to  oj,  Im  CO  >  0,  and  from  cu  to  infinity,  then  we  could  apply 
the  Laplacian  to  the  second  integral  and  the  resultant  integral 
converges  at  Infinity  without  assuming  (^l).   We  could  then  show 
that  those  forms  satisfy  the  various  equations.   We  omit  the 
proofs  as  there  is  nothing  new  in  them.   In  brief,  (57)  and  (58), 
appropriately  interpreted,  are  the  solution  to  our  original  system 
of  equations  without  (^). 

When  f  (v)  and  g(v)  are  spherically  symmetric  P(z,k) 
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and  G(z^k)    are   independent    of  k,    so   fehat 


(59) 


E(r,t)    =    - 


rt     ^ 


vg(v )dv 


r/t 


+ 


Im   r  i   r  ^2   G(z)f/Fr^  sin  (  (r+::t  )/fIFT)    ^   cos  (/FT¥T(r+zt ) )  \ 


^TT^   I J 


2 


r/t 


r/t 


dz  C(z)fj^FrFr  sin((r-zt)XTFT)  ,  cos  (KF(^(r-z_^ 
\  -p  r         > 


In  (59)  in  the  last  two  contour  integrals  the  contours  must  be 
the  same  near  infinity  so  that  the  integral  of  the  sum  of  the 
integrands  may  converge. 

As  is  well  known^   ,  any  spherically  symmetric  distribution 
function  is  stable.  Thus  we  may  use  the  methods  of  Section  IV  to 
obtain  asymptotic  expansions  of  E(r,t).   The  expansion  in  the 
limit  r/t  =  v  fixed,  t  — ^  oo^is  very  complicated  and  none  too 
illuminating,  but  with  little  effort  one  may  show  that 


(60) 


E(r,t) 


h(v) 


In  the    limit   t   — ^  oo,    r   fixed,    we   obtain 


(61) 


E(r,t 


aa)rr        cos(cDt 
o   P  ^    P 


67r/l? 


(V^g) 


3/2 


-+    0(l/t2) 


Ki 


VII.   Two  Dimensional  Green's  Function 

This  section  duplicates  much  of  the  material  of  the  last 
section  for  the  two  dimensional  case.   From  transform  techniques 
we  are  led  to  consider  the  expressions 


(62)   E(r,t) 


Re 
2.2  J 


an  k 


dz  G(z,k)yF(z,k)'sln()/F(z,k)(k-r-zt)) 
k-  r 


+ 


27r^t 


Re 


dn  k  G(z,k) 


k-r/t 


and 


(6:5')  f(r,v,t)  - 


,bf    (v) 
^   Re    ° 


27r  m 


r  r 

an  k 


Sv 


dz  G(z,k) 

k-  r 
t 


cos(/F(z,k)'(k-r-zt)) 
(k-v-z) 


+  g(v)6(r-vt). 

Now  all  vectors  have  two  components  and  dfi  Is  the  "surface 
element"  on  the  "two-sphere",  I.e.,  dfi  =  d0 .   In  the  process  of 
verifying  one  of  the  differential  equations  we  obtain  another 
form  for  the  unusual  Integral  appearing  In  (62).   Since  the 
Integrand  In  (63')  Is  not  properly  defined^ we  consider  an  alter- 
nate form  which  is  meaningful: 
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(63) 


;  (r,v,t 


fSf  (v)   ^ 
Re !  — •  kdQ 


277  m  J  hv 


J 


dz  log(k-v- 
k-  r 


X 


dz  ^ 


,k)cos()/F(^,k)'(k-r-zt)) 


2iY   m 


Re    df)  k 


^f  (v) 


o 


■  k-r 


k-  T' 


a7 


G(^,  t^)log(k-v  -   ^ 


+  g(v ) 6 (r-vt ) . 

By  the  exact  parallels  of  the  arguments  presented  previously 
we  might  easily  show  that  the  Green's  function  is  well-defined, 
that  it  satisfies  the  initial  condition  (2a)  and  the  boundary 
condition  (2b).   Clearly,  E(r,t)  is  curl  free,  so  that  (ic) 
holds,  and  after  the  usual  transformation  of  variables,  r'  =  r-vt, 
and  t'  =  t,  equation  (la)  also  follows. 

If  we  set 


(64)      ^(r,t)  = 


1 


Re 


^   /  k  •  r   ^  \ 
df^  k  G(^,  k)  , 


27r  t    J 


then  we  derive  easily  that 


;65)      V-E(?",t)  -  h 


ire 


f(F,7,t)(d7)=  V  £(r.t)+^g(^) 

t 


From  the  modified  form  of  the  representation  of  G(z,k),  we  learn 


(66) 


r 


£(F,t)  =^Re  \  an   kfk.  Mz^j  iog(k.(7  -  p(d7) 


TTt 


av 
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With  the  help  of  standard  integrals  we  obtain 


Re 


^  A        A  A  A 

dfikk   log   k-a   =   iraa  +   cll    , 


where   II    is   the   unit   dyad   and   c    is    some   constant.      Hence    (66) 
becomes 


(67) 


eir,t 


2e  1= 
or  .. 


(dv)g(v)log(v    -   ^) 


From  (67)  and  (65)  we  see  that  (lb)  is  satisfied. 

Since  there  does  not  seem  to  be  a  simple  form  for  the 
electric  field  when  g(-v)  and  f  (v)  are  functions  of  |v|  only, 
it  is  difficult  to  obtain  both  of  the  asymptotic  expansions 

obtained  previously.   However,  in  the  cylindrically  symmetric 

(4) 
case  --  for  which  the  distribution  function  is  again  stable^   - 

it  is  possible  to  obtain  an  asymptotic  expansion  in  the  limit 

r  fixed,  t  — ^  <».   in  this  case  the  one  dimensional  methods  may 

A 

be  applied  to  the  contour  integral  and  then  the  k  integral  may 
be  performed.   Following  this  technique  we  obtain 


(68) 


E(r,t) 


^         Sm  GD  t 

2   a        p 
-TTOD  a  r  r  ^^^- 


1 


+  0(-p) 
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Footnotes 

(1)  For  various  approaches  to  this  problem  see,  for  Instance, 

L.  Landau,  J.  Phys .  (USSR)  10,  25  (19^0), 
N.  G.  Van  Kampen,  Physlca  21_,  9^9  (1955), 
and  the  comprehensive  study  of  Backus, 

G.  Backus,  J.  Math.  Phys.  _1,  I78  (196O). 

(2)  A  preliminary  account  of  this  work  was  presented  at  the 
November  I961  meeting  of  the  Division  of  Plasma  Physics 

of  the  American  Physical  Society  at  Colorado  Springs,  Colorado, 

(3)  See  Courant-Hilbert ,  Methods  of  Mathematical  Physics,  Vol.  II 
Interscience-Wiley  (1962),  to  appear  shortly. 

(4)  See,  for  instance,  G.  Backus,  op.  cit. 
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